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Extinction transition of bacteria under forced rotation is analyzed in pie geometry. Under 
convection, separation of the radial and the azimuthal degrees of freedom is not possible, and the 
linearized evolution operator is diagonalized numerically. Some characteristics scales are compared 
with the results of recent experiments, and the "integrable" limit of the theory at narrow growth 
region is analyzed. 
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The time evolution of bacterial colonies on a petri-dish 
has been studied recently both theoretically and exper- 
imentally The colony is relatively simple biologi- 
cal system, and its basic component, an individual bac- 
terium, involves only "elementary" biological processes 
like diffusion, food consumption, multiplication, death 
and perhaps some interaction like chemotaxis. Studies 
with some bacteria strains have reported a wide vari- 
ety of complex pattern formation, in most cases due to 
competition for food resources and chemical interaction. 
With uniform, not-exhaustible background of nutrients 
and without the presence of mutations and chemical sig- 
naling, these simple strains are suppose to invade a region 
of nutrient rich agar in the form of a front propagating 
with some typical velocity, known as the Fisher front 

Biological problems of colony growth in inhomogeneous 
environment and under forced convection has been mod- 
eled recently by Nelson and Shnerb Q and by Dahmen, 
Nelson and Shnerb d). These studies have focused on 
the spectral properties of the linearized evolution oper- 
ator, which becomes non-Hermitian in the presence of 
convection An experiment designed to test these 

predictions, has been carried out recently by Neicu et al 

i- 

In the experiment, a colony of Bacillus subtilis bac- 
teria is forced to migrate in order to "catch up" with a 
shielded region on the the petri-dish, where all the other 
parts of the dish are exposed to an ultra violet (UV) 
light, which (under the experimental conditions) makes 
the unshilded bacteria immotile. It was assumed that the 
adaptation of the bacterial colony to the moving shielded 
region has nothing to do with information processing or 
mutual signaling in the colony, but is attributed solely 
to the combined effect of "dumb" diffusion of individ- 
ual bacteria and the larger growth rate under the shel- 
ter. Theoretically, it was predicted that the adaptation 
of the colony to the moving environment fails as the drift 
is faster than the Fisher front velocity, and in this case 
the colony lags behind the shelter and an extinction tran- 
sition takes place. 

In order to get the essence of the theory, let us con- 
sider a one-dimensional example, where bacteria are dif- 



fusing on a line parametrized by x, and are subject to 
some environmental heterogeneity which implies fluctu- 
ating growth rate. If the bacteria diffuses, multiply and 
are forced to migrate with some convection velocity v, 
the differential equation which describes the evolution of 
the colony is. 



9c _ ^ d'^c 
dt dx^ 



dc 
dx 



+ a{x)c — hc^ 



(1) 



where a{x) is the local growth rate. When the hostile 
environment outside the "oasis" causes the immediate 
death of any bacteria, and inside the oasis there is some 
positive growth rate. 



■i{x) ^ 



a < x < 
oo elsewhere. 
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If there is no drift, the linearized version of this prob- 
lem is equivalent to the (imaginary time) evolution of 
quantum particle in an infinite potential well, and is de- 
termined by the eigenvalues of the evolution operator, 
which gives a colony localized on the oasis if it has some 
minimal width (which scales like the width of the Fisher 
front), xq > TTy^D/a. The introduction of a drift term 
into (|l|) is compensated by the "gauge" of the evolution 
(Liouville) operator eigenfunctions 



(pn — sin{mrx/xQ) e*"^/'^"^ sin(n7ra;/a;o) 
together with the eigenvalues "rigid" shift 



r„(«) = T^{v = 0) 
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The theory, thus, predicts an extinction transition as all 
the eigenvalues of (^ becomes negative, i.e., as Vc = 
2y/ aD — 0{1 / xq) , which is the Fisher velocity ||l^. Right 
above the extinction transition, only the largest growth 
eigenvalue (the "ground state") is positive, hence the 
nonlinear interaction between eigenmodes [the term —bc^ 
in (0)] are suppressed at the transition, and the analysis 
is focused on the ground state of the linearized operator. 

In the experiment Q , part of a petri-dish was shielded 
from the UV source, and then this shield was given a 
constant angular velocity with respect to the petri-dish. 
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The corresponding convection velocity v{r) = ujr was 
chosen to interpolate between zero (at the rotation axis) 
and about 2vc at the edge of the dish. It turns out that 
the colony indeed fails to keep rotating with the shield 
at about half the radius. On the other hand, the velocity 
profile for the bacterial density c(r, oj) does not equili- 
brate on the time scales of the experiment (~ 3 days) . In 
this paper, I consider the differences between the one 
dimensional system (|l|,^) and the actual experimental 
setup. In particular, the two-dimensional nature of the 
experiment and the effect of radial diffusion are consid- 
ered explicitly. 

In order to capture the essential physics using the sim- 
plest geometry, the same extinction transition is consid- 
ered on a pie, i.e., a section of the two dimensional disc [ 
See Fig. (1)]. Although the shielded region of the exper- 
iment Ql was not in that shape, it turns out that even in 
this simple geometry there is a coupling between the ra- 
dial and the azimuthal degrees of freedom, and the spec- 
trum becomes "chaotic" when convection takes place. 
The results for this case are, accordingly, relevant also 
to the more complicated geometry of the experiment. 

The basic equation for the bacterial growth problem 
on a non-uniform substrate, in the absence of mutation 
and chemical interactions, is [pi: 



9c(x, t) 
dt 



DV^cix, t) + a(x)c(x, t) + V • Vc - bc^ 



(5) 

With no convection and homogeneous, positive a this 
equation supports Fisher front propagation with veloc- 
ity 1\jDa. The experimental situation corresponds to 
D ^ 10~^cm^/s and a ~ 10~^/s. The Fisher velocity is 
of order 0.1 — Ifim/s, as has been found experimentally. 
The Fisher width, which is the characteristic scale of spa- 
tial correlations in their system, is \/D/a ~ 10~^cto, 
much smaller than the petri-dish radius of few centime- 
ters. 

In cylindrical geometry, Eq. (|^) takes the form, 
dc{r, e, t) 



dt 



= DV-'cir, 9, t) + a{e)c{r, e,t)+vVc~- bc^ 



(6) 




FIG. 1. Pie geometry: the growth rate in the shaded 
region is a, while any bacteria outside this area are dying in- 
stantly due to the UV light. The shaded region is then rotated 
at angular velocity u). 



and for rotating petri-dish the convection term is , 

(7) 

Pie geometry is defined by. 



aiO)^ ^ '^'f'\ (8) 
— oo elsewhere 



i.e., we have absorbing boundary conditions |1J]: 

c(r,6lo,t) = c(r,0,t) = 0. (9) 

As for the boundary conditions on the petri-dish edge at 
r = _R, it is reasonable to take Von-Neumann boundary 
and to impose the no-slip condition on the bacterial den- 
sity at the surface. However, the data |Q| seems to indi- 
cate extinction of the colony at the edge of the dish. This 
is perhaps due to the fact that the width of the bound- 
ary layer (which is expected due to the no slip condition) 
is about the Fisher width, which has been shown above 
to be very small. Accordingly, we further simplify the 
problem by using. 



c{R,0t) = 0. 



(10) 



Dropping the term — bc^ at Eq. (||), one has the linearized 
evolution operator, and for the no-drift (w = 0) case, the 
density of bacteria at time t is given by: 



c{r,et) = ^ A™,„e'^-^'"-0„,„(r,0), 

m,n 

with the eigenstates of the evolution operator, 



(11) 
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4'ni.n{r,0) ^ri,n.nJ^{r/\ D/rm,n)sm{——), (12) 

V On 



where rjm^n is normalization factor 
2 1 



i?V^j[^+i](i?//D7rW)^ 



(13) 



and the constants v4„j_„ are determined by the initial den- 
sity distribution c(r, 6,t — 0). The eigenvalues of the Her- 
mitian problem are: 



= D 



R 



(14) 



where j^^j^^ is the m-th zero of the corresponding Bessel 
function. 
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Let us get an order of magnitude estimate for the time 
scales which are relevant to the experiment Q . The char- 
acteristic times needed for the "ground state" to control 
the system is given by the typical difference between two 
eigenvalues. In our case, since the first zeroes of the 
Bessel functions are of order 1, the times involved are 
~ For an experimental system with R ~ 0.01m 

and D ^ 10~^°m^/s, the typical relaxation times are 
O(lO^sec) ^ 11 days, which is larger than the typical 
time of the actual experiment. 

Now I look at the non-Hermitian case, where w ^ 0. 
Unlike its Cartesian analogy , no simple gauge solves 
the problem, and separation of variables is impossible. 
Spanning the space of normalizable functions by a set of 
Hermitian eigenstates, the perturbative term udg mixes 
both quantum numbers m and n. The matrix elements 
of the convection term are: 



(n, mjcj • delk, I) = 2aj_R^7„ 



where •jnmki is: 



k + n = even 
k + n = odd 



2fen 
.^•2 _^-2 ' ln,m 



Vk,lln 



(15) 



(16) 



and 



Inmki^ J2^{jZ,/eoy) J^ijkr^/eoy)ydy- (17) 

Jo ° 

In order to get the eigenvalues and eigenfunctions at fi- 
nite angular velocity one should diagonalize the full non- 
Hermitian Liouville operator, and the extinction transi- 
tion takes place as the ground state (smallest) eigenvalue, 
Fi^i, becomes larger than the growth rate a on the pie. 

As the rotating system is not integrable, it should be 
studied numerically using some computer diagonaliza- 
tion of the linearized evolution operator. Essentially, one 
should look at the ground state of this operator, since 
this state dominates the system close to the extinction 
transition. 

This numerical analysis, however, may lead to erro- 
neous results if the continuum limit is not taken care- 
fully. In the most general case, a discretized version 
of a model with local growth rate and hopping between 
sites may be realized numerically as a matrix, where the 
growth rates are the coefficients on the diagonal and the 
hopping process gives the off-diagonal terms. As any 
hopping term is positive semi-definite, the only negative 
terms are the local growth rate, and for any finite matrix, 
by adding an appropriate multiplication of the unit ma- 
trix, one may get a positive semi-definite matrix with the 
same eigenvectors. Perron-Frobenius theorem then 
implies that ground state should be a nodeless, positive 
eigenvector. There is a simple physical interpretation to 
this result: since the ground state dominates the sys- 
tem at long times, and the number of bacteria should 



not become negative, Perron-Frobenius theorem should 
hold. Diagonalizing numerically the evolution operator, 
one may get a ground state with nodes, which is physi- 
cally impossible. 
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FIG. 2. Contour Plot of the bacterial density at 6o — n. 
(upper panel) '-j^jj^ ~ 0, (middle); j^j^ = 10 and (lower 
panel) at -jjj^ ~ 30. 
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AO small enough. If the effective discretization is given 
by A6 = Oq/ti, one should truncate the matrix (16) only 
as 



D ' 



(18) 



Although (|l^) seems to indicate that the numerical di- 
agonalization of ( [l5| ) becomes simpler as 6o — > 0, this, 
in fact, is not the case. As the eigenvalues of the un- 
perturbed problem are related to the zeroes of the cor- 
responding Bessel functions, and the rotation operator 
admits matrix elements only between eigenvalues related 
to Bessel functions of different order, it is much simpler to 
diagonalize (|l5|) as >> 0. As Bq ^ 0, the higher m ze- 
roes of any Bessel of order n are smaller than the m = 1 
zero of the n + 1 state and the condition ( |T^ ) implies 
the diagonalization of an infinite matrix. Accordingly, I 
present here the numerical results for the case = i"- 
This situation does not coincide with the experimental 
conditions at Q , but there seems to be no prevention to 
perform the same experiment with large shielded area. 

In Fig. (2), contour plots of the ground state for dif- 
ferent angular velocities are shown. One may identify 
clearly the large deviations from the ground state from 
its shape a.t ui ~ 0. The largest 100 spectral points for 
any case are shown in Fig. (3). 

Fig. (4) presents the ground state eigenvalue, Fq, in 
units of D/R^, as a function of the angular velocity of 
the dish. The extinction transition takes place as this 
eigenvalue is larger than the growth rate on the "pie" , 
jjj'j^i , as has been found above. 
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FIG. 3. First 100 spectral points [I^ij^j]^) vs. 

Re{^jjjj)] at 00 = TT. (upper panel) -^j^jj = 0, (middle): 
~ 10 and (lower panel) at = 30. 

In order to solve this problem one should carefully take 
the discrete limit of the continuum theory. For our case, 
as is discretizcd to quanta of A0, the hopping rate due 
to diffusion becomes Djir^Affl) and the hopping rate 
due to the drift is iw/Afl. In order to avoid the (physi- 
cally impossible) negative hopping rates, one should keep 
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FIG. 4. Highest eigenvalue, -jjy^, as a function of the an- 
gular velocity, -jjj^, for 6^0 = vr. The extinction happens as 
this eigenvalue is larger than a, the growth rate inside the pie. 
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FIG. 5. Narrow shell geometry, where the two-dimensional 
problem converges to the "integrable" case. 

Let us show now how to get a problem equivalent to 
(|l|,||) on a rotating petri-dish. In order to do that, the 
geometry should be taken on a narrow shell as in Fig. 
(5), i.e., the boundary conditions are, 



c(r, 6*0, =c(r,0,i) = 0. 
c(i?i,0o,t) =c(ii'2,0,i) =0. 



(19) 



with AR = R2 — Ri- In the limit i?i ^ 00 at constant 
n, the asymptotic expansion of the Bessel functions Ji, 
and Yi, at large argument gives the eigenfunctions of the 
unperturbed Liouville operator. 



<Pm nir, 0) « ^ — ■ — sm — — 

' ARVR^ ^ AR 



. . , mra . 
an)sm(— — ), 



(20) 



where the phase ensures the boundary conditions at 

2 2 

i?i and the eigenvalues, Tn^m — D ^^yi , are independent 
of n. The matrix elements of the operator udg are given 
by 



(n, m\uj ■ delk, I) = ujSm,i Ink 



with 



Ink 



k + n = even 
k + n ~ odd — 3- 





2kn 



(21) 



(22) 



where the approximation 

r-fl2 „ 




an)sm(-^ +ak)Vrdr 



. ,mnr , . , Inr 
sm(— — + a„) sm(^^ + a^) dr 



AR 



AR 



(23) 



for << 1 has been used. Accordingly, for any m 
sector, both the diagonal and the off diagonal matrix el- 
ements are identical with the corresponding one dimen- 
sional problem, and the results should be the same. 



In conclusion, the mathematical problem which cor- 
responds to the experiment by |4| has been found to 
be non-integrable, and no simple gauge transformation 
connects the eigenvectors of the static and the dynamic 
problems. The actual critical velocity and ground state 
properties have to be studied numerically, and the limit 
of very narrow sector (^o ~* 0) involves diverging numer- 
ical load. The time scales needed for the ground state 
to dominate the system are larger than the duration of 
the actual experiment, and this explains the observed in- 
equilebration. 

I wish to thank A. KudroUi, D. R. Nelson and K. Dah- 
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